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. Abstract. We state, discuss, provide evidence for, and prove in special cases the 

conjecture that the probabihty that a random tihng by rhombi of a hexagon with side 
lengths 2n+a, 2n+b, 2n+c, 2n+a, 2n+b, 2n+c contains the (horizontal) rhombus with 



coordinates {2n + x,2n + y) is equal to i + ga,b,c,x,y{n){^^)^ ^ (3^), where ga,b,c,x,y{n) 

is a rational function in n. Several specific instances of this "1/3-phenomenon" are 
made explicit. 

(N 
> 

o 

O I 1. Introduction and statement of the conjecture 

Let a, b and c be positive integers, and consider a hexagon with side lengths a, b, c, a, 
Q ! b,c whose angles are 120° (see Figure l.a for an example). The subject of our interest 
is the enumeration of tilings of this hexagon by rhombi (cf. Figure l.b; here, and in the 
sequel, by a rhombus we always mean a rhombus with side lengths 1 and angles of 60° 
and 120°). 

As is well-known, the total number of rhombus tilings of a hexagon with side lengths 
a, b, c, a, b, c equals 

a b c . .J 



i=i j=i k=i 



(This follows from MacMahon's enumeration |n. Sec. 429, g — >■ 1; proof in Sec. 494] of 



all plane partitions contained in an a x 6 x c box, as these are in bijection with rhombus 
tilings of a hexagon with side lengths a, b, c, a, b, c, as explained e.g. in [Q.) 

The problem that we are going to address in this paper is the problem of enumerating 
rhombus tilings of a hexagon which contain a given fixed rhombus. Since the total 
number of rhombus tilings of a given hexagon is known, thanks to MacMahon's formula, 
we may ask equivalently the question of what the probability is that a rhombus tiling 
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a. A hexagon with sides a, b, c, a, b, c, b. A rhombus tiling of a hexagon 

where a = 3, 6 = 4, c = 5 with sides a, b, c, a, b, c 

Figure 1. 

of a hexagon that is chosen uniformly at random (to be precise, it is the tiling which 
is chosen at random, while the hexagon is given) contains a given fixed rhombus. (For 
example, we may ask what the probability is that a randomly chosen rhombus tiling 
of the hexagon with side lengths 3, 5, 4, 3, 5, 4, shown in Figure 4, contains the shaded 
rhombus. At this point the thick lines are without relevance.) 




Figure 2. 
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If this question is asked for an "infinite" hexagon, i.e., if we imagine the 2- dimensional 
plane being covered by a triangular grid (each triangle being an equilateral triangle; see 
Figure 2; at this point shades in the figure should be ignored), and ask the question of 
what the probability is that a particular rhombus formed out of two adjacent triangles 
(for example the shaded rhombus in Figure 2) is contained in a randomly chosen rhom- 
bus tiling of the plane (that is compatible with the triangular grid, of course), then 
there is a simple argument which shows that this probability is 1/3: Let us concentrate 
on one of the two adjacent triangles out of which our fixed rhombus is formed. (In 
Figure 3 we have enlarged the chosen rhombus. It is composed out of the triangles 
labelled and 1. We are going to concentrate on the triangle labelled 0.) This triangle 
is adjacent to exactly three other triangles. (In Figure 3 these are the triangles labelled 
1, 2 and 3.) In a rhombus tiling this triangle must be combined with exactly one of 
these to form a rhombus in the tiling. Hence, the probability that a random tiling will 
combine the triangle with the particular one to obtain the fixed rhombus is 1/3. 




Figure 3. 

For a (finite) hexagon however, we must expect a very different behaviour, resulting 
from the boundary of the hexagon. The probability that a particular rhombus is con- 
tained in a random tiling will heavily depend on where the rhombus is located in the 
hexagon. (This is for example reflected in the asymptotic result of Cohn, Larsen and 
Propp 0, Theorem 1].) In particular, we must expect that the probability will usually 
be different from 1/3. 

Rather surprisingly, JohnQ [0, bottom of p. 198] and Propp [|T^ |T8|, Problem 1] made 



the empirical observation that in a hexagon with side lengths 2n — 1, 2n — 1, 2n, 2n — 1, 
2n — 1, 2n the probability that the central rhombus is contained in a random tiling 
is exactly 1/3, the same being apparently true in a hexagon with side lengths 2n, 2n, 
2n — 1, 2n, 2n, 2n — 1. These facts were proved by Ciucu and the author Cor. 3] 
and, independently, by Helfgott and Gessel |]^, Theorem 17]. In fact, more generally, 

^In fact, in pl| the problem of finding the probabihty that, given a hexagonal graph, a chosen fixed 
edge is contained in a randomly chosen perfect matching of the graph is dealt with. The motivation 
to consider this problem is that such hexagonal graphs serve as models for benzenoid hydrocarbon 
molecules. The above probability is called Pauling's bond order. It measures how stable a carbon- 
carbon bond (corresponding to the fixed edge) in a benzenoid hydrocarbon molecule is. 

It is well-known that this problem is equivalent to our tiling problem. The link is a bijection between 
rhombus tilings of a fixed subregion of the infinite triangular grid (such as our hexagons) and perfect 
matchings of the hexagonal graph which is, roughly speaking, the dual graph of the subregion (see e.g. 
[p^ ; "roughly speaking" refers to the little detail that the vertex corresponding to the outer face is 
ignored in the dual graph construction). 
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in both papers the probabihty that in a hexagon with side lengths A^, A^, M, N, N, 
M, N ^ M mod 2, the central rhombus is contained in a random tiling is expressed 
in terms of a single sum, from which the 1/3- result follows on simplification of the 
sum. These results were generalized in two directions. On the one hand, Fulmek and 
the author found a single sum expression for this probability for any rhombus on 
the (horizontal) symmetry axis of the hexagon. On the other hand, Fischer gave 
a single sum expression for the probability that the central rhombus is contained in a 
random tiling of a hexagon with arbitrary side lengths (i.e., with side lengths a, b, c, 
a, b, c). Some further single sum expressions for probabilities of "near-central" rhombi 
to be contained in a random tiling of a hexagon with sides A^, A^, M, A^, A^, M have 
been derived in Theorem 2] and [0. Finally, in complete generality, Fischer P, 
Lemma 2] and Johansson [jlO|, (4.37)] found triple sum expressions for the probability 
that a fixed (not necessarily central or near-central) rhombus is contained in a random 
tiling of a hexagon with side lengths a, b, c, a, b, c. (These two triple sum expressions 
are completely different from each other.) 

The purpose of this paper is to report a curious manifestation of the fact that "in 
the limit" the probability that a particular rhombus is contained in a random tiling is 
1/3. Roughly speaking, it seems that the probability equals 

1/3 plus a "nice" expression. 

To make this precise, we need to introduce a convention of how to describe the position 
of a rhombus in a given hexagon. First of all, without loss of generality, we may restrict 
our considerations to the case where the fixed rhombus is a horizontal rhombus (by 
which we mean a rhombus such as the shaded ones in Figures 2-4), which we shall 
do for the rest of the paper. (The other two types of rhombi are then covered via a 
rotation by 120°, respectively by 240°.) In order to describe the position of a rhombus 
in the hexagon, we introduce, following the following oblique angled coordinate 
system: Its origin is located in one of the two vertices where the sides of lengths b and c 
meet, and the axes are induced by those two sides (see Figure 4). The units are chosen 
such that the grid points of the triangular grid are exactly the integer points in this 
coordinate system. (That is to say, the two triangles in Figure 4 with vertices in the 
origin form the unit 'square.') Thus, in this coordinate system, the bottom-most point 
of the shaded hexagon in Figure 4 has coordinates (5,4). 

With this convention, we have the following conjecture. It extends an (ex) conjecture 



by Propp [|17[ |l^. Problem 4]. 



Conjecture. Let a, b, c, x and y be arbitrary integers. Then the probability that a 
randomly chosen rhombus tiling of a hexagon with side lengths 2n + a, 2n + b, 2n + c, 
2n + a, 2n + b, 2n + c contains the (horizontal) rhombus with bottom-most vertex 
{2n + x,2n + y) (in the oblique angled coordinate system) is equal to 



1 „ . . /2nV / /Qn + 2 



3 + J /^3^^J forn>no, (1.1) 

for a suitable uq which depends on a, b, c, x and y, where fa,b,c,x,y{n) is a rational 
function in 



^This statement is clearly equivalent to the statement in the abstract. The form (1.1) of the 
expression is more convenient in the subsequent listing of special cases. 




Figure 4. The oblique angled coordinate system 



As is shown in Section 2, for any specific a, b, c, x, y the corresponding formula 
for fa,b,c,x,y{n) can be worked out completely automatically by the use of a computer 
(given that the Conjecture is true, of course). We have in fact produced a huge list 
of such formulas, of which we list a few selected instances below. As we explain in 
Section 3, any of these is (at least) a "near-theorem," in the sense that it could be 
proved automatically by the available multisum algorithms, provided there is enough 
computer memory available (and, thus, will at least be a theorem in the near future). 
Also in Section 3, we elaborate more precisely on which of these are just conjectural, 
and which of them are already theorems^. However, we do not know how to prove the 
Conjecture in general, that is, for generic values of a, b, c, x, and y (cf. Section 3 for a 
possible approach). 

Here is the announced excerpt from our list of special instances of the Conjecture: 

/-1,-1.0,-1,-iH = /2,2,i,2,i(ri) = forn > 1, (1.2) 

f2,l,l,2,lin) = /2,l,l,l,l(n) = /l,2,l,2,l(^) = /l,2,l,l,o(^) 

= /-i,o,o,o-i(^) = /-i,o,o -1 -i(^) = /o -i,o,o,o(^) = /o -1,0-1 -i(^) = forn > 1, 

(1.3) 



^For the convenience of the reader, we have marked conjectures by an asterisk in the equation 
number. 
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= /l,l,0,l -l(^) = fl,l,0,l,l{^) 
= /l,l,l,0-l(^) = /l,l,l,0,lH = /l,l,l,l,o(^) = /o,2,0,l,oH 

= /2,0,0,0,l(^) = /2,0,0,l,o(^) = /2,0,0,l,l(^) = /2,0,0,2,o(^) = ^ for 71 > 1, (1.4) 

/i,iAO -iH = /i,i,o,o,oH = /2,o,i,i,o(n) = ^ forn > 1, (1.5*) 

2 

/o,2,i,i,o(ri) = -- for n > 1, (1.6) 

/i,MAoH = for n > 1, (1.7*) 
4 

/i,i,o,i,o(^) = 3 for n > 1, (1.8) 
4 

/4,3,l,3,2(n) = /4,3,l,4,2(n) = - for 71 > 1, (1.9*) 

/oAOAoW = -^^^ for ^ > 1, (1-10) 

/o,o,i,o,o(n) = forn>l, (l.U) 

, , 2(2n+l)(3n + 2)Un + 5) 

= 3(n+ l)^(6n + 5) " ^ ^' '"^^ 

, , "in^ + ISn^ + 12n + 1 , , 
/..o,3.-.(n) ^ 6(.+ 1)^(2.-1) f""^^- 

6[n + l)''(rz + 2)(6rz + 5) 
„ , , (37z + 2)(2n2 + 4n + l) , 

/-i,5A.-iH = 3(7.+ l)> + 2) " ^ ^' (^-^^ ) 

•/ 10,3,0,1,4 6(n+l)^(n+2)2(rt+3)2(n+4)^(n+5)(2n-3)(2n-l)(6n+5)(6n+7)(6n+ll) 

x(176n9+3080n8+21692n7+74546n6+102578n5-73279n4-362598n3-283977n2+24762n+55440) 

for 71 > 2. (1.16*) 
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2. How are these conjectures and results discovered? 



Point of departure for all tliese discoveries is an observation by Propp |1^, Prob- 
lem 4]: He conjectured that the probability that a randomly chosen rhombus tiling of a 
hexagon with all side lengths equal to contains the "near-central" rhombus (this is 
the rhombus with bottom-most vertex (A^, A^) in the oblique angled coordinate system) 
is equal to 1/3 plus a "nice" formula in Should this observation be true, then the 
Mathematica program Rate[| ("Rate!" is German for "Guess!"), respectively its Maple 
equivalent GUESS0, will find the formula, given enough initial terms of the sequence. 

Let us see how this works in the case that is odd. For the generation of the 
probabilities, Propp used the programs vaxmapleQ and vaxmacsQ, which are based on 
the evaluation of determinants of large (if though sparse) matrices. However, since 
then triple sum formulas have been found by Fischer Lemma 2] and Johansson flO 



(4.37)], which allow to generate these probabilities much more efficiently. We choose to 
use Fischer's formula. We state it below. 

Theorem. Let a, b and c be positive integers, and let (x, y) be an integer point such 
that < X < b+a — 1 and 1 < y < c+a — 1. Then the probability that a randomly chosen 
rhombus tiling of a hexagon with side lengths a, b, c, a, b, c contains the (horizontal) 
rhombus with bottom-most vertex {x,y) (in the oblique angled coordinate system) is 
equal to 



j — l\ fc + i + x — y — 2\ /b + s — x + y — 1 
: — 1 J \ b + s — X — 1 

(6+ l),_i (c+ l)i_i {b + c + i)j- 



+ l)r \S —\) \ X — 1 J \ b + S — X — 1 

' i=\ j=l s=l 



(z-l)!(6 + c+l),_i 

We now program this formula in Mathematica. 
Mathematica 2.2 for DOS 387 
Copyright 1988-93 Wolfram Research, Inc. 
In[l] : = F [a_ ,b_,c_,x_,y_] :=c! /Pochhammer [b+1 , c] * 
Sum [Sum [Sum [ (-1)^ (i+s) *Binomial [j-1 , s-1] * 

Binomial [c+i+x-y-2 ,x-l] *Binomial [b+s-x+y-1 ,b+s-x-l] * 
Pochhammer [b+1 , s-1] *Pochhammer [c+1 , i-1] * 

Pochhammer [b+c+i , j-i] / ( j-i) ! / (i-1) ! /Pochhammer [b+c+1 , s-1] , 



(2.1) 



■^Commonly, by a "nice" formula one means an expression which is built by forming products and 
quotients of factorials. A strong indication that one encounters a sequence {aN)N>o for which a "nice" 
formula exists is that the prime factors in the prime factorization of do not grow rapidly as N 
becomes larger. (In fact, they should gro w linearly.' 



^written by the author; available from http : //radon, mat .univie . ac. at/People/kratt; see |13 



Appendix A] for an explanation of how the program works. 

^written by Frangois Bcraud and Bruno Gauthier; available from 



tittp : //www-igm . univ-mlv . fr/ "gauthier . 

''written by Greg Kupcrbcrg, Jim Propp and David Wilson; available at 



tittp : / /math .wise . edu/~propp/sof tware .html 



written by David Wilson; also available at ittp: //math. wisc.edu/~propp/sof tware. html. 
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{s,l,j}] ,{j,l,a}] ,{i,l,a}] 

Now we generate the first eleven values of these probabilities for N = 2n + 1 and 
subtract 1/3 from them. 

In[2] := Table[F[2n+l,2n+l,2n+l,2n+l,2n+l]-l/3,{n,l,ll}] 

4 3 2000 245 296352 142296 43188288 

Out [2]= { , , , , , , , 

105 143 138567 22287 33393355 19126225 6743906935 

759169125 15365378600 55469016746 805693639296 

) } 

135054066707 3067656658059 12280863528759 195909013434965 
Next we load Rate, and apply Rate's function Ratekurz to the sequence of numbers. 
In[3] := ( (rate.m 
In[4] := Apply [Ratekurz, 7o2] 

2 

(1 + 2 il) (2 + 3 il) (4 + 3 il) 

4 Product [ ,{il, 1, -1 + iO}] 

2 

(1 + il) (5 + 6 il) (7 + 6 il) 

Out [4]= { } 

105 

The program outputs a formula which generates the terms of the sequence that was 
given as an input. The formula is written as a function in iO, i.e., we must replace iO 
by n. In more compact terms, the formula can be rewritten as 



1 /2n\ ^ / /6n + 2 



3\nJ / \3n + 1 



(2.2) 



(It should be observed that this expression is exactly the one which features in ( |1 . 1| ) . ) 

At this point, this formula is of course just a conjecture. It has however been proved 
in 0, Corollary 7, (1.9)]. 

Being adventurous, one tries the same thing for other choices of the parameters a, 
b, c, X and y. Very quickly one discovers, that a similar phenomenon seems to occur 
for any choice 2n + a, 2n + b, 2n + c for the side lengths and (2n + x,2n + y) for the 
coordinates of the bottom-most point of the fixed rhombus, where a, b, c, x, and y are 
fixed integers. Although the (conjectural) expressions that one finds need not be "nice" 
anymore in the strict sense above, it is at worst polynomial factors in n that appear 
in addition. Moreover, one also realizes soon that division of such an expression by 
the expression in (|2.2| ) apparently always results in a rational function in n, i.e., the 
Conjecture in Section 1 is discovered. 

Let us see just one such example. We choose a hexagon with side lengths 2r;, + 2, 
2n+l, 2n, 2n + 2, 2n-|-l, 2n, and (2n + 3,2n — 1) for the coordinates of the bottom-most 
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point of the fixed rhombus. Then we obtain the following numbers for n = 1, 2, 
The reader should note that we immediately divide the expression ( |2.2| ). 

In [5] : = Table [ (F2 [2n+2 , 2n+l , 2n , 2n+3 , 2n-l] - 1/3) / (Binomial [2n , n] ^3/ 

Binomial [6n+2,3n+l] ) ,{n,l,15}] 

35 43 307 593 337 1585 2339 1099 4483 5921 
Out [5]= {-(--), ~, , , , , , , , , 

12 54 480 1050 648 3234 4992 2430 10200 13794 

2545 9649 11987 4891 17731 
) , , , , } 

6048 23322 29400 12150 44544 

By having a brief glance at this sequence, it seems that the first term is "alien," so let 
us better drop it. 

In [6] := Drop[%,l] 

43 307 593 337 1585 2339 1099 4483 5921 

Out [6]= {--, , , , , , , , , 

54 480 1050 648 3234 4992 2430 10200 13794 

2545 9649 11987 4891 17731 
) , , , , } 

6048 23322 29400 12150 44544 

By the discussion above, this should be a sequence which is given by a rational function 
in n. Therefore is suffices to apply Rate's Rateint (which does just rational interpola- 
tion, in contrast to Ratekurz, which tries several other things, and which is therefore 
slower). 

In [7] := Apply [Rateint, 7,] 

2 3 

35 + 60 iO + 30 iO +4 iO 

Out [7]= { } 

2 

6 (2 + iO) (1+2 iO) 

Again, the program outputs the formula as a function in iO. Since initially we dropped 
the first term of the sequence, we must now replace iO by n — 1. 

In[8]:= Factor [7,/. iO-)n-l] 

2 3 

l+12n+18n +4n 

Out [8]= { } 

2 

6 (1 + n) (-1 + 2 n) 
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Hence, if the Conjecture in Section 1 is true, /2,i,o,3 must be the expression given 
in the output Out [8] . (Thus, we have discovered Eq. (1.13).) Again, at this point, this 
is just a conjecture. 

The Equations ( p..2D -(1.16) in Section 1 are all found in the same way. 



3. Discussion: How to prove the conjecture? 

A possible approach to prove the Conjecture in Section 1 is to start with the ex- 
pression (|2.1|) (or with the alternative expression [1^, (4.37)]), replace a hj 2n + a, b 
by 2n + b, c by 2n + c, x by 2n + x, y by 2n + y, and by some manipulation (for ex- 
ample, by applying hypergeometric transformation and summation formulas) convert 
it into the form ( |2.2|) . Everybody who has some experience with manipulating bino- 
mial/hypergeometric sums will immediately realize that this is a formidable task. In 
particular, it seems a bit mysterious how one should be able to isolate "1/3" from the 
"rest." In any case, I do not know how to prove the Conjecture in this manner, nor in 
any other way. 

On the other hand, as we explained in Section 2, for any specific values of a, b, c, x, 
and y, it is routine to find a conjectural expression for the rational function fa,b,c,x,y{n) 
(given that the Conjecture is true). In turn, once such an expression is available, it 
can (at least in principle) be verified completely automatically. For, what one has to 
prove is the equality of the expressions (p.l|), with the above replacements, and (|1.1|), 
where fa,b,c,x,y{n) is the explicit rational function found by the computer. That is to 
say, one has to prove that a certain triple sum equals a closed form expression. Clearly, 
this can be done (again, at least in principle) by the available multisum algorithm^, by 
using the algorithm to find a recurrence in n for the expression ( ^.Ij ) , and subsequently 
checking that the expression (|1.1|) , with the computer guess for fa,b,c,x,y{n), satisfies 



the same recurrence. Unfortunately, in any case that I tried, the computer ran out of 
memory. 

However, as we already mentioned in the Introduction, in some cases formulas in form 
of single hypergeometric sums are available. If one is in such a case then one would 
proceed as in the above paragraph, but one would replace the multisum algorithm by 
Zeilberger's algorithm]^ (see |1T5| , 16, 23, 24|). The advantage is that, in contrast to 



the multisum algorithm, Zeilberger's algorithm is very efficient. At any rate, in any 
case that I looked at in connection with our problem, the Zeilberger algorithm was 
successful. That is to say, if I am allowed to somewhat overstate it, whenever one is 
in a case where a single sum formula is available, one has a theorem (i.e., Zeilberger's 
algorithm will prove that the empirical found rational function fa,b,c,x,y{n) does indeed 
satisfy the Conjecture for all values of n). 



^The first (theoretical) algorithm for proving multisum identities automatically was 
given by Wilf and Zeilberger |2^. A considerable enhancement and speedup was ac- 
complished by Wegschaider who combined the ideas of Wilf and Zeilberger with 
ideas of Verbaeten pO|. Wegschaider's Mathematica implementation is available from 



tittp : //www .rise. uni-linz . ac . at /research/ combinat/risc/ software. 

^"^A Maple implementation written by Doron Zeilberger is avail- 



able from ittp : //www. math, temple . edu/~zeilberg ; a Mathematica imple- 



mentation written by Markus Schorn and Peter Paule is available from 



tittp : //www .rise. uni-linz . ae . at /research/ eombinat/rise/ software 
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For the sake of completeness, we list the vectors ( y) for which single sums 

are available for fa,b,c,x,y{n). Clearly, it suffices to restrict c to and 1. (All other values 
can be attained by shifts of n.) 

(A) By [|, Theorem 1, (1.2)]: (2a' + 1, 2b' + 1, 0, a' + b' + 1, a'), for integers a' and b'. 

(B) By [|, Theorem 1, (1.3)]: (2a', 26', 1, a' + b', a'), for integers a' and 6'. 

(C) By [I, Theorem 2, (1.4)]: (2a' + 1, 26' + 1, 1, a' + b' + 1, a' + 1) and (2a' + 1, 26' + 
1, 1, a' + 6' + 1, a'), for integers a' and 6'. 

(D) By p. Theorem 2, (1.5)]: (2a', 26', 0, a' + 6', a') and (2a', 26', 0, a' + 6', a' - 1), for 
integers a' and 6'. 



(E) By [|, Theorem 1] 

(F) By [|, Theorem 2] 

(G) By Theorem 3] 



(a', a', 0, 2x' + 1, x'), for integers a' and x'. 
(a', a', 1, 2x', x' — 1), for integers a' and x'. 
(2a', 2a', 1, 2a', a' + 1) and (2a', 2a', 1, 2a', a' — 1), for an integer 



a'. 



(H) By 0, Theorem 4]: (2a' + 1, 2a' + 1, 0, 2a' + 1, a' + 1) and (2a' + 1, 2a' + 1, 0, 2a' + 

1, a' — 1), for an integer a'. 
(I) By ]3, Theorem 5]: (2a', 2a', 0, 2a', a' + 1) and (2a', 2a', 0, 2a', a' - 2), for an integer 



a'. 



(J) By 0, Theorem 6]: (2a' + 1, 2a' + 1, 1, 2a' + 1, a' + 2) and (2a' + 1, 2a' + 1, 1, 2a' + 
1, a' — 1), for an integer a'. 

Thus, choosing a' = 6' = in (A), we see for example that the expression for 
/o,o, 1,0,0 ('^) given in (|1.11|) is in fact a theorem. For, by Theorem 1, (1.3) in [^] with 
a = b = 2n, c = 2n + 1, the probability that a randomly chosen rhombus tiling of 
a hexagon with side lengths 2n, 2n, 2n + 1, 2n, 2n, 2n + 1 contains the (horizontal) 
rhombus with bottom-most vertex (2n, 2n) can be written in the form 



SUM(n):^g^;i(^HH2- 
^ ' ^ (2n + l)4„ \nj\nj 



n + 3/2)fc (2n + V 



k 



■{n + k + 2)„_,_i {2n + k + 2)„_fc_i /^^^^" \\ (3.1) 

[n — k — ly. 

Next we take it as an input for Zeilberger's algorithm (we are using Zeilberger's Maple 
implementation here): 

|\'^/| Maple V Release 4 (Uni Wien) 

._|\| l/l_. Copyright (c) 1981-1996 by Waterloo Maple Inc. All rights 

\ MAPLE / reserved. Maple and Maple V are registered trademarks of 

( ) Waterloo Maple Inc. 

I Type ? for help. 

) read ekhad: 

) ezra(zeillim) ; 

z e i 1 1 im ( SUMM AND ,k,n,N, alpha, beta) 
Similar to zeil(SUMMAND,k,n,N) but outputs a recurrence for 
the sum of SUMMAND from k=alpha to k=n-beta . 
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Outputs the recurrence operator, certificate and right hand side. 
For example, "zeillim(binomial(n,k) ,k,n,N,0,l) ; " gives output of 

N-2, k/(k-n-l) ,1 
which means that SUM(n) :=2^n-l satisfies the recurrence 
(N-2)SUM(n)=l, as certified by R(n,k) :=k/(k-n-l) 

) zeillim(2*n*(2*n+l) ! /rf (2*n+l,4*n)*binomial(2*n,n)*binomial(3*n,n)* 
) 2^ (2*n-2) *rf (n+3/2 , k) *rf (2*n+l , k) *rf (n+k+2 , n-k- 1 ) * 
) rf (2*n+k+2,n-k-l)*rf (l/2,n-k-l)/(n-k-l) ! ,k,n,N,0,l) ; 

-1 + N, 1/6 (1 - 2 n + 2 k) (-288 n - 432 n k - 912 n - 1440 n k 

42 32 3 3 22 2 

- 216 n k - 414 n k - 1746 n k - 1096 n - 189 n k - 612 n 

2 3 2 3 2 

+ 36 n k - 907 n k + 48 n k - 152 n - 163 n k + 43 n k - 12 + 4 k 

2 3 / 

+ 32 k + 16 k ) / ((-n + k) n (6 n + 1) (6 n + 5) (2 n + k + 2) 
/ 

2 

(1 + 3 n + 2 n )), 

2 3 n 3 2 
GAMMA (3 n) GAMMA (n + 1/2) 64 (36 n + 60 n + 29 n + 3) 
1/2 

2 3/2 3 2 

(n + 1) Pi GAMMA(n) (6 n + 5) (6 n + 1) n GAMMA(6 n) 



It tells us that the expression SUM(n) in satisfies the recurrence 

(3 + 29n + 60n2 + 36n^) {3n - 1)!^ (2n)! 



2n2 (n + 1)^ {6n + 1) (6n + 5) (n - 1)!^ {6n - 1)! n!^ 

(3.2) 



SUM(n + 1) - SUM(n) 



(The first term in the output, -1 + N, encodes the form of the left-hand side of 
the third term gives the right-hand side. The middle term is the so-called certificate 
which provides a proof of the recurrence.) So it just remains to check that the expression 
( [1.1|) with /o,o,i,o,o(^) as in ( |1.11| ) satisfies the same recurrence and agrees with (|3.1| ) for 
n = 1, which is of course a routine task. 

On the other hand, the expression for /2,i,i,2,i('^) given in (|1.3|) cannot be established 
in the same way by appealing to a special case of one of (A)-(H). Still, it is also 
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a theorem, thanks to the following simple observation: suppose that we consider a 
hexagon with side lengths a, b, c, a, b, c, where a = b, and a rhombus on the horizontal 
symmetry axis of the hexagon. Let us imagine that this rhombus were the one in 
Figure 3 (consisting of the triangles labelled and 1). Let us denote the probability 
that a randomly chosen tiling contains this rhombus by p. Since the rhombus is on 
the symmetry axis, the probability that a randomly chosen tiling contains the rhombus 
consisting of the triangles labelled and 2 is equal to the probability that it contains the 
rhombus consisting of the triangles labelled and 3. Let us denote this probability by 
q. Any tiling must contain exactly one of these three rhombi, hence we have p + q + q = 
p + 2q = 1. Therefore, whenever there is a single sum formula available for p, there is 
also one for q. To come back to our example, the rhombus whose bottom-most point 
has coordinates {2n + 2, 2n + 1) in a hexagon with side lengths 2n + 2, 2n + 1, 2n + 1, 
2n + 2, 2n + 1, 2n + 1, can be seen as such a rhombus consisting of triangles labelled 
and 3, where the bottom-most point of the rhombus consisting of the triangles labelled 
and 1 has coordinates {2n + 1, 2n + 1). This puts us in Case (E), with a' = x' = — 1, 
(to see this one has to replace n by n — 1 in the above coordinatization) , and thus the 
claimed expression for /2,i,i,2,i('^) can be proved in the same manner as we proved the 
expression for /o,o,i,o,o('^) above. 

Again, for the sake of completeness, we list the additional vectors (a, 6, c, x, y) for 
which single sums are available for fa,b,c,x,y{n) by the above observation. 

(C) By 0, Theorem 1]: (2a', 0, 0, a' + 1, a'), (2a', 0, 0, a', a'), (2a', 0, 0, a', a' - 1), (2a', 0, 

0, a' - 1, a' - 1), (0, 2a', 0, a' + 1, 0), (0, 2a', 0, a', -1), (0, 2a', 0, a', 0), (0, 2a', 0, a' - 

1, -1), for an integer a'. 

(D') By [g Theorem 2]: (2a' + 1, 1, 1, a' + 2, a' + 1), (2a' + 1, 1, 1, a' + 1, a' + 1), (2a' + 
1, 1, 1, a' + 1, a'), (2a' + 1, 1, 1, a', a'), (1, 2a' + 1, 1, a' + 2, 1), (1, 2a' + 1, 1, a' + 1, 0), 
(1, 2a' + 1, 1, a' + 1, 1), (1, 2a' + 1, 1, a', 0), for an integer a'. 

(E') By [|, Theorems 1 and 2]: (2a', 0, 0, a' + x', a' - x'), (2a', 0, 0, a' + x', a' - x' - 1), 
(2a' + l, 0, 0, a' + x', a'-x'), (2a' + l, 0, 0, a' + x', a'-x' + l), (0, 2a', 0, a' + x', 2x'-l), 
(0, 2a', 0, a' + x', 2x'), (0, 2a' + 1, 0, a' + x', 2x' - 2), (0, 2a' + 1, 0, a' + x', 2x' - 1), 
for integers a' and x'. 

(F') By [|, Theorems 1 and 2]: (2a', 1, 1, a' + x', a' - x' + 1), (2a', 1, 1, a' + x', a' - x'), 
(2a'+l, 1, 1, a'+x', a'-x'+2), (2a'+l, 1, 1, a'+x', a'-x'+l), (1, 2a', 1, a'+x', 2x'-l), 
(1, 2a', 1, a' + x', 2x'), (1, 2a' + 1, 1, a' + x', 2x' - 2), (1, 2a' + 1, 1, a' + x', 2x' - 1), 
for integers a' and x'. 

(C) By [g Theorem 3]: (2a' + 1, 0, 0, a' + 2, a' + 2), (2a' + 1, 0, 0, a' + 1, a' + 2), (2a' + 

l,0,0,a',a'), (2a' + l,0,0,a'-l,a'), (0, 2a' + l, 0, a' + 2, 0), (0, 2a' + 1, 0, a' + 1, -1), 

(0, 2a' + 1, 0, a', 0), (0, 2a' + 1, 0, a' - 1, -1), for an integer a'. 
(H') By 0, Theorem 4]: (2a', 1, 1, a' + 2, a' + 2), (2a', 1, 1, a' + 1, a' + 2), (2a', 1, 1, a', a'), 

(2a', 1, 1, a' - 1, a'), (1, 2a', 1, a' + 2,1), (1, 2a', 1, a' + 1,0), (1, 2a', 1, a', 1), (1, 2a', 1, 

a' — 1, 0), for an integer a'. 
(F) By §, Theorem 5]: (2a', 0, 0, a'+2, a'+l), (2a', 0, 0, a'+l, a'+l), (2a', 0, 0, a'-l, a'), 

(2a', 0,0, a' - 2, a'), (0, 2a', 0, a' + 2, 0), (0, 2a', 0, a' + 1, -1), (0, 2a', 0, a' - 1, 0), 

(0, 2a', 0, a' — 2, —1), for an integer a'. 
(J') By §, Theorem 6]: (2a' + 1, 1, 1, a' + 3, a' + 2), (2a' + 1, 1, 1, a' + 2, a' + 2), (2a' + 

l,l,l,a',a'-l), (2a' + l, 1,1, a'-l, a'-l), (1, 2a' + 1, 1, a' + 3, 1), (1, 2a' + 1, 1, a' + 

2, 0), (1, 2a' + 1, 1, a', 1), (1, 2a' + 1, 1, a' - 1, 0), for an integer a'. 
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a, b, c, X, y 


1 • / A \ /TT\ 

entry m (Aj-(Hj, 
respectively (C )-(H ) 


—1, —1, 0, —1, —1 


(A), a = b = -1 


2, 2, 1, 2, 1 


(B), 0=0=1 


2, 1, 1, 2, 1 


(F ), a = X = 1 


2, 1, 1, 1, 1 


(F'), o' = 1, x' = 


1, 2, 1, 2, 1 


/-r^l\ f ft 

(F ), a = X = \ 


1, 2, 1, 1, 


(F'), o' = 1, x' = 


-1, 0, 0, 0, -1 


(E'j, a = -1, x' = 1 


—1, 0, 0, —1, —1 


(E ), a = —1, x—{) 


0, -1, 0, 0, 


(E ), o = —1, X = 1 


0, —1, 0, —1, —1 


(E'), a' = —1, x' = 


1, 1, 1, 1, 1 


(C), 0=6=0 


1, 1, 0, 1, —1 


/TT\ / T / 1^ 

(H), 0=6=0 


1, 1, 0, 1, 1 


/TT\ / T / 1^ 

(H), 0=6=0 


1, 1, 1, 0, —1 


(F), o = 1, X = 


1, 1, 1, 0, 1 


( TIN / f\ 

(J ), 0=0 


1, 1, 1, 1, 


(C), 0=6=0 


0, 2, 0, 1, 


(E'), o' = 1, x' = 


2, 0, 0, 0, 1 


/TIN / -1 

(F), o = 1 


z, U, U, i, U 


(^Hj j, o — i, X — U 


2, 0, 0, 1, 1 


(E'), o' = 1, x' = 


2, 0, 0, 2, 


(E'), o' = x' = 1 


0, 2, 1, 1, 


(B), o' = 0, 6' = 1 


1, 1, 0, 1, 


(A), o' = 6' = 


0, 0, 0, 0, 


(D), o' = 6' = 


0, 0, 1, 0, 


(B), o' = 6' = 


3, 3, 0, 3, 1 


(A), o' = 6' = 1 



Table 1. 



Table 1 lists the special cases that we considered in ( |1.2| )- ([L^ ), ( |1.6| ), ( |1.8| ), (|1.10|) - 
( [1.12| ), for which proofs are available, together with an indication from which of the 
Cases (A)-(H), respectively (C')-(H'), these proofs come from. In particular, the ex- 
pression for /o,o,i,o,o(^) in ( |1-11D (together with ( |1.1| )) provides the formula for the values 
of Pauling's bond order in Tabelle 2 of |Tl[| for higher naphtalenes N{p) of odd order 
p = 2n — 1, and the expression for /i,i,o,i,o('^) in (together with ( |1.1| )) provides the 
formula for the values of Pauling's bond order in Tabelle 3 of ||TT| for higher pyrenes 



P{p) of odd order p = 2n — 1. The values in Tabelle 1 of |TI|] for higher benzenes B{p) 
are expressed by (^TTj) with a = b = c = x = y = and a = 6 = c = x = ?/=l, re- 
spectively, with the expression for /o,o,o,o,o('^) given in ( |1.10|) and the one for /i,i,i,i,i(^) 



given in (|1.4|) . (The latter formulas have already been stated in Cor. 7, (1.8) and 
(1.9)]. In fact. Corollary 7 of contains some more evaluations of this kind.) 

Coming back to the original goal, a proof of the Conjecture for arbitrary a, b, c, x and 
y, it may seem that it should be at least possible to achieve this in the Cases (A)-(H) 
and (C')-(H'), where single sum formulas are available. For, for each specific choice of 
a, b, c, X and y out of one of these cases, an identity of the form "single sum = closed 
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form" has to be proved. So one would try to follow the strategy that was suggested in 
complete generality at the beginning of this section: apply some manipulations (using 
hypergeometric transformation and summation formulas, for example) until the desired 
expression is obtained. This task is much less daunting here, since we are dealing now 
with a single sum, not with a triple sum. Moreover, as it turns out, the sums that occur 
are very familiar objects in hypergeometric theory (we refer the reader to |^, |TP|, ^ for 
information on this theory), they turn out to be balanced 4F3-series, respectively very- 
well-poised yFg-series. (For example, the series in ( |3.1| ) is a balanced 4F3-series.) For 
these series there are a lot of summation and transformation formulas known. However, 
and this is somehow mysterious, I was not able to establish any of the theorems that I 
presented here in this classical manner (i.e., without the use of Zeilberger's algorithm), 
not to mention a general theorem for an infinite family of parameters. As already said 
at the beginning of this section, the biggest stumbling block in such an attempt is the 
question of how one would be able to isolate "1/3" from the "rest." So, potentially, 
there is a hierarchy of interesting hypergeometric identities lurking behind the scene 
which has not yet been discovered. 
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